Particle Systems with Stochastic Passing 
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We study a system of particles moving on a line in the same direction. Passing is allowed and 
when a fast particle overtakes a slow particle, it acquires a new velocity drawn from a distribution 
Po(v), while the slow particle remains unaffected. We show that the system reaches a steady state 
if Po(v) vanishes at its lower cutoff; otherwise, the system evolves indefinitely. 
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We are interested in behavior of a system of interacting 
particles moving on the real line in one direction, say 
to the right. The system is endowed with the following 
simple dynamics: (i) Particles move freely between "col- 
lisions" ; (ii) After a collision, or "passing" event, the ve- 
locity of the slow particle remains the same, i> s i ow =const, 
while the fast particle instantaneously acquires some new 
velocity, v ncw > v s \ ow , drawn from the intrinsic velocity 
distribution Po(v). We would like to answer the following 
basic questions about the behavior of the system: Does 
the velocity distribution P(v, t) reach a steady state or 
the system continues to evolve indefinitely? How does 
the average velocity depend on time? etc. 

Our motivation is primarily conceptual, as we want to 
understand non-equilibrium infinite particle systems with 
two-body interactions. Thus we have chosen the simplest 
dynamics - interactions occur only upon colliding, and 
only one particle is affected. The appealing simplicity of 
the model suggests that it might show itself in different 
natural phenomena, and indeed, we originally arrived to 
this model in an attempt to mimic traffic on one-lane 
roads. Somewhat related dynamics were already used in 
modeling voting systems [QB , force fluctuations in bead 
packs , asset exchange processes Q] , combinatorial pro- 
cesses || , continuous asymmetric exclusion processes |1 , 
granular gases 0, and aggregation- fragmentation pro- 
cesses ||. 

Let us first consider discrete velocity distributions. 
Specifically, we assume that both initial velocities and 
new velocities are drawn from the same intrinsic distribu- 
tion Po(v) — ^PjS(v — Vj). For the binary distribution, 
the system does not evolve at all, so the first non-trivial 
case is the ternary intrinsic distribution when the system 
contains slow, moderate, and fast particles. Initially, 

Po(v) = p\5{v - vi) +p 2 S(v - v 2 ) + p 3 S(v - u 3 ). (1) 



We set 



Pi +P2 +P3 = 1, vi < v 2 < v 3 , 



(2) 



without a loss of generality. When the steady state is 
reached, the velocity distribution remains ternary, 

P eq (w) = PiS(v - vi) + q 2 S(v - v 2 ) + q 3 S(v - v 3 ). (3) 

The density of slow particles does not change, while the 
densities q 2 and q 3 of moderate and fast particles differ 



from the initial values. The final densities are found from 
a simple probabilistic argument based on the requirement 
of stationarity. For moderate particles we get 

p 3 (v 2 - vi)q 2 = p 2 (v 3 - Vi)q 3 . (4) 

The left-hand side of Eq. (g) gives the loss in q 2 which 
happens when a moderate particle overtakes a slow par- 
ticle and becomes a fast particle. The right-hand side 
gives the gain in q 2 which takes place when a fast parti- 
cle overtakes a slow particle and converts into a moderate 
particle. Solving Eq. (Q) together with the normalization 
condition, pi + q2 + 33 = 1, we find 



92 = P2 



P2+P3 
P2 + vp 3 '' 



93 = vp 3 ■ 



P2 +P3 

P2 + vp 3 ' 



(5) 



where v — (v 2 — v\)j(v 3 ~ v{). Since v < 1, we have 
q 2 > p 2 and g3 < p 3 . Thus the density of moderate parti- 
cles increases while the density of fast particles decreases. 
Similarly, one can analyze discrete velocity distributions 
with more than three particle species. In all cases (i) 
the system reaches a steady state; (ii) the average ve- 
locity decreases and eventually reaches some finite value; 
(iii) the density of the most slow particle species remains 
unchanged. 

Now we consider a continuous velocity distributions. 
Let [fmin , v max ] be a support of the intrinsic velocity 
distribution Pq(v). By Galilean transform, we can set 
Vmin = without loss of generality. We consider un- 
bounded distributions, w max = 00, although main results 
equally apply to the cases with finite v max . 

The passing rule asserts that when a fast particle over- 
takes a slow particle moving with a velocity u s iow, the 
assignment of the new velocity v occurs with probability 



Po (vKiow) = Pq(v) 



8{V ~ ^slow) 

lZdv'Po(v'Y 



(6) 



Eq. (|6|) guarantees that v > v s \ ow and that the normal- 
ization requirement, J dv Po (v\v s i ow ) = L is obeyed. 

Now we can write a Boltzmann equation for the veloc- 
ity distribution P(v,t): 



dP(v,t) 
dt 



P(v,t) / dv'(y-v')P(v',t) 



(T) 





00 



+ / dv 2 P (v\v 2 ) dvi{vi-V2)P(vi,t)P{v 2 ,t). 

'0 Jv 2 



1 



The first term on the right-hand side of Eq. (Q) describes 
loss in P(v, t) due to collisions with slow particles: Colli- 
sions occur with rate proportional to velocity difference, 
and the integration limits ensure that only collisions with 
slower particles are taken into account. The second, gain, 
term accounts for the increase of P(v, t) due to a random 
assignment of velocity v after collision. 

We could not solve Eq. (fjj) in the general case of an ar- 
bitrary intrinsic velocity distribution Pq(v). Attempts to 
find a solution even for some particularly simple Pq(v), 
e.g., linear, exponential, or uniform, turned out to be 
fruitless as well. Thus we proceed by employing asymp- 
totic, approximate, and numerical techniques. 

We start by looking at the asymptotic behavior of P(v) 
in the small velocity limit. Let v <C u(t), where u(t) is 
the average velocity, 



(8) 



u(t) = (v) = / dvvP{v,t). 
Jo 

Then Eq. (Q) simplifies to 

dP(v,t) 



in 



= P (v)u(t) / dv 2 P(v 2 ,t) 



P(v,t) / dv'(v-v')P(v',t). 
Jo 



(9) 



To probe the small v behavior, we need to know Po(v) 
at v — > 0. Let us consider a family of intrinsic velocity 
distributions that behave algebraically: 



P (v) ~ Av 1 - 1 when 



0. 



(10) 



Now assume that the system reaches the steady state: 
P(v,t)-^ Peq(v) and u(t) — » u cq . Plugging these and 
Eq. (|l^) into Eq. (g), we find that the velocity distribu- 
tion also behaves algebraically in the small velocity limit, 



P eq (v) ~ (n + l)Au cq v^ 



when 



0. 



(11) 



In other words, the steady state velocity distribution 
scales as i> -1 Po(v). Recalling the normalization require- 
ment, J Q dvP(v) = 1, we see that it is possible only 
when fi > 0. Thus our assumption that the system 
reaches a steady state is certainly wrong when \i < 0. 
In this region, we anticipate that the system will evolve 
indefinitely. Note that both the exponential and uni- 
form intrinsic distributions belong to the borderline case 
of /J, = that separates stationary and evolutionary 
regimes; for them an anomalously slow kinetics is an- 
ticipated. 

To probe the behavior of evolving systems, we assume 
that in the long-time limit there is a very small frac- 
tion of "active" particles that move with velocities v ~ 1 
and the vast majority of "creeping" particles that hardly 
move at all. We ignore collisions between active particles 
since their density is very low. We also ignore collisions 
between creeping particles since their relative velocity is 



very small. This picture suggests that only collisions be- 
tween active and creeping particles matter. Thence, the 
velocity distribution of active particles obeys 



dP(v,t) 
dt 



P {v)u{t)-vP{v,t). 



(12) 



Eq. ( |l2| ) may at best describe the evolution process in 
the long-time limit. However, for the sake of tractability, 
we apply it to the whole time range and use the natural 
initial condition P(v,0) = Pq(v). Eq. ( |i"2| ) is an inhomo- 
geneous linear differential equation which is easily solved 
to give 



P(v,t)=P (v)e 



1 



dt'u(t') 



-vt' 



This solution implies 

P(v,t) ~u(i))j- 1 PoW 



for 



v»- t , 



(13) 



(14) 



which resembles Eq. (pd|). 

To close the solution of Eq. ([l3|), we must determine 
u(t). It is possible to plug ( |l3| ) into the definition of the 
average velocity, Eq. (|8|) , and get an integral equation for 
u(t). In the following we use another approach, which is 
technically simpler. Note that the density of active par- 
ticles, J dvP(v,t), is manifestly conserved by Eq. (|lT 
After integration over velocity, Eq. ( |l3|) becomes 



l = Po(*)+ / dt'u{t')P (t-t'), 



(15) 



where Pq is the Laplace transform of the intrinsic velocity 
distribution, 



Po(t) 



dv Pq(v) e 



(16) 



One can guess the long time behavior of P(v) with- 
out actually solving Eq. ([l5]). Let us assume that the 
average velocity varies slowly with t. Then the integral 
on the right-hand side of Eq. Jl5| ) can be estimated as 
u(t) J* dt' Po(t'), which implies 



u(t) 



dt'P (t') 



(17) 



For an intrinsic velocity distribution with an algebraic be- 
havior ( [To| ) in the small-w limit, we have Po(t) ~ 
for large t. Hence L dt' Po(t') ~ for /i < 0, and it 
follows from Eq. ( |l7| ) that u(t) ~ t^. The above deriva- 
tion is quite careless, though the final result is correct. 
Now we derive this result in a more rigorous way. 

The convolution form of the integral in Eq. (^5|) sug- 
gests to apply the Laplace transform once more. It yields 



o s + v Jo s + v 



(18) 



2 



Here 



u(s) 



dt u(t) i 



(19) 



is the Laplace transform of the average velocity. The 
double Laplace transform of the intrinsic velocity distri- 
bution has been simplified: 



dti 



dv Pq (v) e 



dv 

lo Jo Jo s + v 

Thus the Laplace transform of the average velocity is 

-l 



Po(v) 



(20) 



u(s) = -1 + 



s + v 



(21) 



Generally, one cannot obtain more explicit results. Given 
that the above approach describes only the long-time 
asymptotics, let us focus on this regime. To probe the 
long-time behavior, one should determine the small s 
asymptotics of u(s). For algebraic intrinsic velocity dis- 
tributions (|To|) , the asymptotics of (|2(]) reads 



dv 
o s 



Po(v) 



As 11 [ dw — — = ^ — - 

J w + 1 sm(— iTfj,) 



(22) 



This applies for — 1 < fi < (the lower bound comes from 
the normalization requirement, J dv P (v) = 1). Plug- 
ging (H) into @ yields 



u(s) 



sin(-TT^) 
An 



for 



0. 



(23) 



and by making the inverse Laplace transform, we finally 
arrive at 



u(t) 



sin(— 7T ju) 
AnT{l + n) 



t" for t -> oo. 



(24) 



This result agrees with the asymptotics we naively de- 
rived earlier. 

A special consideration is required for the borderline 
case of /i = 0. For concreteness, consider the exponen- 
tial intrinsic distribution, Po(v) = exp(— v). Its double 
Laplace transform reads 



where 



dv — 

o s 



E 1 (a) 



e s E 1 (s), 



dx ■ 



n x 

is the exponential integral. As a result, Eq. (|2l| ) becomes 

1 



u(s) = -1 + 



se s ^i(s) 



(25) 



Using the well-known asymptotics of the the exponential 
integral |§, Ei(s) = - In s - 7 + 0(s), (where 7 = 0.5772 
is Euler's constant), we transform Eq. (pa) into 



u{s) 



1 



O 



1 



s(lns + 7) V ms , 
Performing the inverse Laplace transform gives 

u(t) — » - — for t — > 00. 
mi 



(26) 



(27) 



To summarize, for the family of intrinsic velocity distri- 
bution with algebraic behavior in the small v limit (|ic|), 
our predictions for the long-time asymptotics of the av- 
erage velocity u(t) are: 



{const, for fi > 0; 
(lnt)-\ (oth = 0; 
V 1 , for -1< fi < 0. 



(28) 



To check the validity of asymptotic predictions and, 
more generally, to see if the mean-field theory is appli- 
cable at all, we perform molecular dynamics simulations 
and solve the Boltzmann equation (|7j) numerically To 
sample distinct regimes predicted in (E8) we consider the 
intrinsic velocity distribution 



Po(v) 



T( f i + 1) 



(29) 



with /I = 1,0,-1/2. 

In molecular dynamics simulations, we place N parti- 
cles onto the ring of length L = N so that the average 
density is equal to one. Most of our simulations are per- 
formed for N — 5 • 10 4 particles, but we also simulated 
twice larger system and found no appreciable difference. 
Initially, particle velocities are randomly drawn from the 
distribution Pq{v). The model is updated according to 
the collision-time-list algorithm suggested in Ref. Jl(J. 

To solve the Eq. ([?]) numerically, we use Euler's time 
update with both uniform and non-uniform grid; in the 
latter case, we take vn — (N j N max ) A v max velocity grid 
with v max = 15 and N max — 300 — 500. Integrals on 
the right-hand side of Eq. ([?]) are calculated using the 
trapezoid rule; time increment St = OA was found to be 
suitable for all three Pq(v). 

The main conclusion is that results of molecular dy- 
namics simulations and numerical solutions of the mean- 
field equation are virtually identical (see, e.g., Fig. 1). It 
confirms our assumption that the system remains well- 
stirred and no appreciable spatial correlations develop. 
Additionally, Fig. 1 shows that for Pq(v) = ve~ v , the 
approach of P(v,t) to the steady state is non-uniform 
in velocity. This is caused by the obvious fact that for 
any finite time, the velocity distribution P{y, t) must still 
vanish at the lower cutoff as Po(v) does. In other words, 
the steady state (|ll]) is reached outside the "bound- 
ary layer", v 3> while within the boundary layer, 
v -C v*(t), the velocity distribution continues to evolve. 
The threshold velocity v*(i) is estimated by evaluating 
the first, leading term in the right-hand side of Eq. (0): 



3 



f-ip ~ v^ +1 P, which implies - f-W+i). The width 
of the boundary layer shrinks with time but the boundary 
layer still exists ad infinitum. 

4.0 I 1 1 1 1 1 1 




Fig. 1. Plot of P{v,t) at t = 200 for P (v) = ve~ v : sim- 
ulation result (o), numerical solution ( — ). The dashed line 
shows simulation result for P(v,t) at t = 16000. 



Figs. 2-4 plot the average velocity vs. time for the 
intrinsic velocity distributions (|29|) with /i = 1,0, —1/2, 
respectively. We find good agreement with the theoreti- 
cal prediction of Eq. ( pq ) when \i > 0. On Fig. 5 the plot 
of the local exponent a(t) = dln[ii(t)]/<f ln[t] vs. i~ 1//2 is 
shown for /i = —1/2. The results of extrapolation of a(t) 
to the t — > oo limit are not contradicting the theoretical 
prediction, a — 1/2. 

2.0 | 1 1 1 1 1 1 
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Fig. 2. Plot of u(t) vs. time t for P (v) = v e~ v . 
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Fig. 3. Plot of exp[l/u(t)] vs. time t for P (v) = e~ v . 
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Fig. 4. Log- log plot of u(t) vs. time t for Pa(v) = 
(nv)~ 1 / 2 e~ v : molecular dynamics results (o) and numerical 
solution ( — ). 
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0.00 0.10 0.20 0.30 0.40 0.50 
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Fig. 5. Plot of the local exponent a(t) vs. for 
Po(v) = (t™)- 1 /^"". 

In summary, we have shown that the fate of the sys- 



tern of passing particles is determined by the behavior of 
the intrinsic velocity distribution near its lower cutoff: If 
Pq(v) vanishes in this limit, the system reaches a steady 
state; otherwise, the evolution continues forever. Com- 
parison between solutions of the mean-field Boltzmann 
equation and results of molecular dynamics simulations 
suggests that the mean-field theory description is exact. 
It will be interesting to confirm this result rigorously. 
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